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ABSTRACT 



Aims. We present the first theoretical SN la light curves calculated with the time-dependent version of the general purpose model 
atmosphere code PHOENIX. Our goal is to produce light curves and spectra of hydro models of all types of supernovae. 
Methods. We extend our model atmosphere code PHOENIX to calculate type la supernovae light curves. A simple solver was imple- 
mented which keeps track of energy conservation in the atmosphere during the free expansion phase. 

Results. The correct operation of the new additions to PHOENIX were verified in test calculations. Furthermore, we calculated theo- 
retical light curves and compared them to the observed SN la light curves of SN 1999ee and SN 2002bo. We obtained LTE as well as 
NLTE model light curves. 

Conclusions. We have verified the correct operation of our extension into the time domain. We have calculated the first SN la model 
light curves using PHOENIX in both LTE and NLTE. For future work the infrared model light curves need to be further investigated. 

Key words, stars: supernovae: general - radiative transfer - methods: numerical 



1. Introduction 

All types of supernovae are important for the role that they play 
in understanding stellar evolution, galactic nucleosynthesis, and 
as cosmological probes. Type la supernovae are of particular 
cosmological interest, e.g., be cause the dark energy was discov- 
ered w ith Type la supernovae (Rie ss et al.ll998HPerlmutter et alJ 
Il999h . 

In dark energy studies, the goal now is to characterize the 
nature of the dark energy as a function of redshift. While there 
are other probes that will be used (gravitational lensing, baryon 
acoustic oscillations), a JDEM or Euclid mission will likely con- 
sider supernovae in some form. In planning for future dark en- 
ergy studies both from space and from the ground, it is im- 
portant to know whether the mission will require spectroscopy 
of modest resolution, or whether pure imaging or grism spec- 
troscopy will be adequate. Several purely spectral ind i cators 
of peak luminosity have been proposed (Nugent et al. 1995; 
[ Hachinger et alj|2006t iBongard et alj|2006t iBronder et al.ll2008: 
iFolev et aTl 120081 : IChotard et al.l 120 1 ll) . What is required is an 
empirical and theoret ical comparison of both light curve shape 
luminosity indicators (lPskovskiill977tlPhillipsfl9 93; Riess et al] 
1 1 996t iGoldhaber et al.l 200 1 ) and spectral indicators. 

To make this comparison one needs to know more about the 
physics going on in a supernova explosion and to be able to cal- 
culate light curves and spectra self-consistently. Thus, we need 
to extend our code to time-dependent problems. While our pri- 
mary focus is on Type la supernovae, the time-dependent radia- 
tive transfer code is applicable to all types of supernovae, as well 
as to other objects, e.g., stellar pulsations. 

In the following we present the methods we used to imple- 
ment time dependence. First we focus on solving the energy 
equation (first law) and then present our first theoretical LTE 



light curves of type la supernova events. In a further section, we 
present results of NLTE model light curve calculations. 

2. Energy solver 

In our previous paper dJack et al.ll2009l) . we presented our ap- 
proach to calculate light curves of type la supernovae events. We 
implemented a simple energy solver into our general purpose 
model atmosphere code PHOENIX where we kept track of the 
overall energy change of the radiating fluid and the energy ex- 
change between the matter and radiation. Simple test light curves 
confirmed that our approach worked correctly. We now present 
a new approach to calculate theoretical light curves of SNe la. 
Again, we are using a simple solver where we keep track of the 
energy conservation. In the new approach, we consider only the 
energy density of the material. We assume free expansion and do 
not solve the equations of hydrodynamics. In order to obtain the 
result for the new time step, we use an explicit scheme. We apply 
this scheme and compute time steps until radiative equilibrium 
is reached. These are the points in our light curves. 

The direct change of the energy density of the material 
considering absorption and emission of radiation and energy 
deposition by gamma ra ys is given by equation (96.7) in 
Mihalas&Mihalas(1984) 





de 


d 




p 


df 







[cxE - Anrf) + pe, 



(1) 



where p is the density, p the gas pressure and e is the energy 
density of the material. The quantities of the radiation field are 
X which is the absorption coefficient, r\ is the emission coeffi- 
cient and E — —J is the radiation energy density with the mean 
intensity J. To obtain E, the radiative transfer equation in spheri- 
cal symmetry is solved including special relativity. All additional 
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energy sources are put in e such as the energy input from gamma 
ray deposition. This equation represents the first law of thermo- 
dynamics for the material. The change of the energy density of 
the material depends on the coupling of matter and radiation 
field, the absorption of gamma-radiation and positron annihila- 
tion energy, the expansion work and the change of the ionization 
and excitation of matter. 

Dividing equation 1 by the material density p we obtain 



de _ 1 r 
dt~pj 



(cyE - Ann) dA - p he. 

df p 



(2) 



The radiation energy density is given by E — — J, where J is the 
mean intensity. Using this, we obtain for the radiation term 



(c X E - An^dA = An \ ( X J - j/)(L1 = An { X { J - S ))dA, (3) 



where S 



is the source function. All these quantities can 



be derived from the solution of the radiative transfer equation. 
The term of the change of the energy density by the radiation is 
therefore given by 



Q= XA(JA~S A )dA. 



(A) 



Another change of the energy density is due to the work W 
done by the adiabatic expansion of the SN la atmosphere. For a 
discrete step this work is given by 



W = p[---). 

\P2 Pi , 



(5) 



The expansion is assumed to be homologous. Since we solve the 
energy equation for the matter, we do not include the radiation 
pressure work. Since the system is radiation dominated there is 
the possibility of numerical inaccuracies in coupling the matter 
and radiation only by Q. For the calculation of the new radii 
and densities a s well as a discussion about the accuracy of this 
assumption see I Jack et al.l (120091) . 

Including all energy changing effects, the new energy density 
of the material after a discrete time step Af is explicitly given 
by 



/ 1 l\ An C 
e 2 = e\ - p\ + — Af x(J -S)dA + eAt, 

\P2 Pi) P J 



(6) 



where e\ is the old energy density. 

To obtain all the needed quantities we have to solve the 
spherically symmetric special relativistic radiative transfer equa- 
tion. The advantage of our new approach is that we do not have 
to iterate for each time step. In our previous paper the calcula- 
tion of each time step involved an iteration process to obtain the 
new matter temperature. With the new approach we can calcu- 
late the new temperature for the next time step directly from all 
the known quantities. 

The translational energy density of the material is given by 



em 



2p 



3R r 



3 N A kT 
2 m 



(7) 



with the mean molecular weigh p and the universal gas constant 
R. The gas pressure is represented by p and the density by p. T 
stands for the temperature of the material. This equation of the 
energy density is now used to determine the new temperature 
after the next time step. 



During the first phase of the SN la envelope evolution, the 
material of the atmosphere is hot and, therefore, highly ionized. 
The energy change due to ionization and excitation changes of 
the atoms present in the SN la atmosphere cannot be neglected. 
PHOENIX already solves the equation of state (EOS), where all 
the excitation and ionization stages of the present atoms and 
molecules are included. Using the EOS, we obtain the overall 
energy density of the material by the sum of the ionization en- 
ergy e ions and the translational energy e, mns 

(8) 



-trans ' c tons- 



Hence, the energy density change of the material goes into 
a change of the translational energy and the ionization energy, 
which both depend on the temperature. Therefore, we obtain the 
new temperature by an iteration scheme. The matter density at 
the next point in time is determined by homologous expansion. 
A first temperature guess is used, and the EOS is solved to ob- 
tain the ionization energy density. Combined with the transla- 
tional energy density, the overall energy density is computed. 
This is checked against the target energy density, which we ob- 
tained equation [6] If the obtained energy density is incorrect, a 
new temperature guess is made. This new temperature guess is 
obtained by assuming a linear dependence of the energy density 
and temperature. The current temperature guess is iterated to the 
target energy density. It takes about 5-10 iteration steps to deter- 
mine the new temperature. If the EOS delivers the correct target 
energy density, the new temperature of the next time step has 
been found. The accuracy relative of the energy density in this 
iteration process is set to 10~ 5 . 

2. 1 . y-ray deposition 

The maximum of the light curve of an SN la event is observed 
around 20 days after explosion. Causing this later maximum of 
the light curve of an SN la event is the energy release into the 
envelope caused by the radioactive decay of 56 Ni and its also ra- 
dioactive decay product 56 Co. Therefore, this energy deposition 
has a strong influence on the energy change of the SN la enve- 
lope structure. Hence, the energy deposition due to radioactive 
decay has to be taken into account for the calculation of the SN 
la envelope evolution. 

The energy deposition due to the y-rays emitted by radioac- 
tive isotopes needs to be computed by a radiative transfer solver 
for the y-rays. In this work, we solve the y-ray depositio n with 
the as sumption of a gray atmosphere for the y-rays. Qefferv 
(1998) did a detailed study of the y-ray deposition and pointed 
out that this is an adequate approach to calculate y-ray depo- 
sition in SN la atmospheres. In the decay of a 56 Ni nucleus, a 
y-photon is emitted with an energy of 2.136 MeV. The 56 Co nu- 
cleus decays to an 56 Fe nucleus and emits a y-photon, which 
has an energy of 4.566 MeV. In the decay of 56 Co about 19% 
of the energy is released by positrons. The positrons are as- 
sumed to be locally trapped. They annihilate by emitting two 
photons each with an energy of 512 keV, which has to be taken 
into account for the energy deposition calculation. The opac- 
ity is considered to be constant and a pure abs orption opacity , 
meaning that no scattering is assumed. As in iJeffervl ([1998), 
K y = 0.06 < Z/A > cm 2 g~' was chosen as the opacity. < Z/A > 
is the proton fraction, which is counts both bound and free elec- 
trons, since the electron binding energy is small compared to 
the energy of gamma rays. The energy deposition into the atmo- 
sphere per unit time is given by 

e = An- J, (9) 
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where J is the mean intensity, which has been obtained by solv- 
ing the gray radiative transfer for the y-rays. This obtained en- 
ergy deposition has to be taken into account for the calculation 
of the overall energy change. 

2.2. Adaptive time step procedure 

The timescale for energy changes in the SN la envelope will 
change during the evolution of the light curve. In order to save 
computation time, the light curves have to be calculated with the 
optimal time step size for each phase of light curve evolution. 
Therefore, we implemented an adaptive time step routine to de- 
termine the optimal time step size for the current time step. 

The energy change Ae of the energy of the material e may be 
approximated by 



Ae 



• e = At ■ (Q + e) , 



(10) 



where Q is the energy change of the interaction with the radia- 
tion, and e is the energy deposition by the gamma rays. The en- 
ergy change due to the expansion is ignored in this case. On the 
one hand, this energy change depends on the new matter density 
after the time step, which is unknown because it depends on the 
size of the time step itself. Furthermore, the energy change be- 
cause of the expansion is small compared to the changes caused 
by the energy transport and the energy deposition by y-rays. The 
idea of the adaptive time step procedure is to limit the energy 
change to a prescribed amount of the energy of the material. 
Thus, rewriting equation [10] we obtain the time step size Af for 
the current time step by 



Af = 



■ X, 



(11) 



where x is the introduced limiting energy change factor. The fac- 
tor x ranges between x, m „ and x max , which mark the largest and 
smallest allowed energy change. These are input parameters for 
the adaptive time step procedure. The time step size is calculated 
for every layer, and the minimum time step size of all layers is 
used for the energy solver. 

Each time the adaptive time step procedure is called, it 
checks if the energy changes of the time step before were too 
big or could have been bigger. If the minimum time step size is 
in a different layer, the same value of x is kept for the next time 
step. If the minimum is in the same layer and the sign of the en- 
ergy change does not change, the previous time step might have 
been too small. Thus, the factor x is increased for the following 
time step. If the sign changes, the last time step might have been 
too large, therefore, the factor x is decreased. This means that 
for each time step, the allowed energy change is adapted and the 
factor x is updated to get the optimal time step during the whole 
evolution of the SN la atmosphere. 



3. Test calculations 

All new implemented physical processes of the simple energy 
solver have to be tested. For the test atmosphere, the at mosphere 
struct ure and abundances of the W7 deflagration model (Nomoto 
1 19841) are used. The atmosphere structure is expanded to a point 
in time of 10 days after the explosion. The densities and radii are 
determined by free homologous expansion and can be computed 
easily. To perform the test calculations, we obtained an initial 
temperatur e structure with the PH OENIX temperature correction 
procedure (Hauschi ldt et al.l l2003). This temperature structure is 



5.0x10 4 



hydro solver 
temperature correction 



1x1 4 2x10 4 3x1 4 

expansion velocity [km/s] 

Fig. 1. Temperature structures of a test atmosphere obtained with 
the energy solver and the PHOENIX temperature correction pro- 
cedure. 



the result of a simple approach that does not represent the tem- 
perature structure of an SN la precisely. However, for the test 
calculations for our energy solver, this simple temperature struc- 
ture is sufficient. With this initial atmosphere structure, the en- 
ergy solver is applied for different test cases. All contributions 
to the energy change are tested separately. 

3.1. Energy transport 

In this section, the energy transport through the atmosphere 
is tested. The energy solver considers only the energy change 
caused by emission and absorption of radiation, where the result 
of the radiative transfer equation is needed. All other influences 
are neglected. As a first test, we check how the initial tempera- 
ture structure changes if the energy solver is changing the SN la 
envelope. As the initial atmosphere structure is already in radia- 
tive equilibrium, the energy solver should not change the tem- 
perature structure significantly, because it also pushes the atmo- 
sphere towards a radiative equilibrium state. 

In Fig. Q] a comparison of the temperature structure of the 
energy solver to the result of the temperature correction proce- 
dure is shown. The differences in the temperature structure for 
most layers are less than 1%. But the temperature differences of 
the inner layers are clearly higher. These differences arise in the 
PHOENIX temperature correction result which produces a spike 
in the temperature structure. This is likely due to the bound- 
ary condition in the temperature correction. Hence, the resulting 
temperature structure obtained with the energy solver is more 
accurate. Here, the temperature structure is smooth. In order to 
obtain an atmosphere in radiative equilibrium, the energy trans- 
port part of the energy solver can be used instead of the tempera- 
ture correction procedure. The main problem is that about a few 
hundred time steps are needed to obtain the resulting atmosphere 
structure in radiative equilibrium, while the temperature correc- 
tion needs fewer iteration steps and is, therefore, significantly 
faster. 



3.2. Expansion 

In a further test calculation, the expansion part of the energy 
solver is checked. The only energy change considered is the adi- 
abatic cooling due to free expansion of the SN la envelope. The 
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Fig. 2. Light curve of a test atmosphere that is just expanding Fig. 4. Light curve of a test atmosphere with energy deposition, 
and adiabatically cooling. Because of the additional energy, the luminosity increases with 

time. 
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Fig. 3. Initial and final temperature structure of a test atmosphere 

which is just expanding and adiabatically cooling. Fi 8- 5 - Model H g ht curve of a realistic test scenario. 



80 



energy deposition by y-rays or an energy change due to energy 
transport is disabled. For this test case, the expectation is that 
the atmosphere should just cool down, so the temperature of the 
atmosphere and the observed luminosity should be decreasing. 

The observed luminosity is shown in Fig. [2] The observed 
luminosity of the SN la atmosphere decreases in time. The tem- 
perature structure of the first and the last time step is plotted in 
Fig. [3] The adiabatic expansion has cooled the atmosphere ev- 
erywhere, and the new temperature structure is now significantly 
lower than the initial one. 

3.3. Energy deposition 

To test the energy deposition by the radioactive decay of nickel 
and cobalt into the SN la atmosphere, a test case is considered, 
where only this y-ray deposition is calculated with the energy 
solver. The energy change due to free expansion and energy 
transport are neglected to see the direct effect of the additional 
energy put into the test SN la atmosphere. 

The results of the energy deposition by y-rays calculation 
with the energy solver is shown in Fig. |4] where the light curve 



is plotted. Due to the energy added to the atmosphere, the lumi- 
nosity seen by an observer increases with time. 

3.4. Realistic test scenario 

After all single effects have been tested, we now calculate a test 
case where all effects are considered for the solution of the en- 
ergy solver. Again, the same initial temperature structure is used 
and we start our calculation at day 10 after the explosion. Free 
expansion as well as energy deposition and energy transport are 
active for this computation. 

The observed light curve is shown in Fig. [5] The luminosity 
increases because of the energy input from deposition of energy 
from the y-rays produced by decay of nickel and cobalt. It takes 
a certain time, until the whole atmosphere has relaxed to this 
new condition. The atmosphere is then in radiative equilibrium 
state, and the luminosity stays constant. The initial and final tem- 
perature structure are plotted in Fig. |6] The energy input caused 
by radioactive decay has increased the temperature of the whole 
atmosphere. The atmosphere is heated by y-ray deposition in 
the inner part of the atmosphere. Due to this additional energy, 
the luminosity of these layers increases and the heat is radiated 
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Fig. 6. The initial temperature structure for the static solution of 
the equation of radiative equilibrium and the relaxed temperature 
structure obtained 1 minute later by advancing the energy solver 
in the test scenario. All influences on the SN la atmosphere are 
considered. 



away and absorbed by the surrounding layers. This energy trans- 
port takes care that the deposited energy is moving through the 
whole atmosphere so that the temperature increases everywhere 
and the additional energy from the radioactive decay is radiated 
away towards the observer. The atmosphere is then in radiative 
equilibrium. 



4. SN la model light curves 

We now present our first theoretical SN la light curves ob- 
tained with our extensions to the general purpose model atmo- 
sphere code PHOENIX and compare them to observed SN la light 
curves. The online supernov a spectrum archive (SUSPECT) 
dRichardson et al.ll200lL 20020 provides numerous of observa- 
tions of different types of supernovae. For this work, the ob- 
served light curves of SN 2002bo and SN 1999ee are used to 
compare them to our results of model light curves. Photometric 
light curv e observations of S N 2002bo in different photomet- 
ric bands dBenetti et al.ll2004 have been obta ined. SN 1999ee 
also has observed spec tra (lHamuv et al.ll2002l) and photometry 
dStritzinger et al.ll2002l) . 

4.1. Method 

The energy solver is now applied to calculate synthetic light 
curves of SNe la. The SN la light curve evolution is calculated 
during the free expansion phase. For the initial model atmo- 
sphere structure, the results of the explosion calculation of other 
groups are used as the input structure. Each layer has a certain 
expansion velocity, which does not change during the evolution, 
because homologous expansion is assumed. We start the model 
light curve calculation a few days after explosion. In the first few 
days the SN la envelope is optically thick and compact. Another 
point is, in the first few days the SNe la light curves are quite 
faint, and there are almost no observations of this early time 
phase that have been obtained. Thus, for the model light curve 
calculation, it is adequate to start the light curve calculation a 
few days after the explosion. The initial structure is given by 



http://suspect.nhn.ou.edu/ 



the result of the explosion model simulation. The results of the 
explosion model give the expansion velocities, density structure 
and the non-homogeneous abundances of all chemical elements 
present in the SN la envelope. The envelope expands for the first 
few days by assuming homologous expansion. The radii are de- 
termined by the expansion velocities and time after explosion. 
For the first temperature structure guess, the PHOENIX temper- 
ature correction procedure is used to obtain a initial tempera- 
ture structure, which is in radiative equilibrium. Alternatively 
we could calculate RE models starting at day one. 



For the computation of an SN la model light curve, we let 
the energy solver change the obtained initial model atmosphere 
structure. The atmosphere structure adapts to the new conditions 
caused by y-ray deposition and other energy effects. After a cer- 
tain time, the atmosphere eventually reaches the radiative equi- 
librium state. A typical time step of a energy change in the model 
atmosphere is about 0.1s. It takes about 500 time steps to reach 
radiative equilibrium depending on in which evolution phase the 
SN la is. For the later phase after the maximum of the SN la 
light curve, fewer time steps are needed. We use this first ra- 
diative equilibrium atmosphere structure to calculate a more de- 
tailed spectrum, where more wavelength points are used. This 
model spectrum is then used to obtain the first point of the model 
light curve for each band by using the filter functions described 
below. It certainly would require too much computation time 
to calculate a whole light curve evolution by using the typical 
time steps of about 0.1s. Therefore, to obtain the next point of 
the light curve, big time steps are computed. For these big time 
steps, the atmosphere is only expanding. This means that nei- 
ther energy deposition by y-rays or energy transport through the 
atmosphere is considered for the solution of the energy solver. 
After half a day computed with big time steps, the whole en- 
ergy solver changes the atmosphere structure again and the next 
point in the light curve is obtained after the atmosphere structure 
moves back to radiative equilibrium. Even with the use of these 
big time steps, the radiative transfer for a light curve of 50 days 
has to be solved about 10,000 times. For normal model atmo- 
sphere calculation in PHOENIX with the temperature correction 
procedure, the radiative transfer equation has to be solved about 
100 times, at most. Therefore, the computation of a whole SN 
la model light curve is very computationally expensive, because 
the solution of the radiative transfer equation has to be obtained 
too many times. It takes a huge amount of computation time to 
calculate even simple SN la model light curves, even when using 
only a few wavelength points. Basically, the SN la model light 
curve is a curve consisting of half day points, where the atmo- 
sphere is in radiative equilibrium. During the later phase after 
maximum light, the big time steps have been performed for one 
or even two days, as the energy changes is the SN la atmosphere 
become smaller. 



As described above, a model spectrum is calculated with the 
obtained atmosphere structure, which is in radiative equilibrium 
at each half day point of the light curve. To obtain a model light 
curve for different photometric bands, filter functions were used 
to calculate the luminosity in different bands. In Fig. [7] an SN 
la model spectrum and the filter functions for the U, B, V, R 
and I bands a re sho wn. These filter functions are described in 
lHamuv et alJ ([1992). With these filter functions, SNe la model 
light curves can now be obtained for these five different photo- 
metric bands. 
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plot of filters 




wavelength [A] expansion velocity [km/s] 



Fig. 7. A typical model spectrum of an SN la together with the 
filter functions for the U, B, V, R and I band. These are used to 
determine the luminosity in each band for the model light curves. 



4.2. Light curves of LTE models 

For the first calculations of theoretical light curves, the model at- 
mosphere of the SN la is considered to be in LTE. For a first ap- 
proach to obtain model light curves this is adequate. Another rea- 
son for the assumption of LTE is that model light curve calcula- 
tions with an atmosphere treated in NLTE use more computation 
time. In order to obtain an SN la model light curve in reasonable 
computation time, a model atmosphere in LTE is a necessary 
assumption. For the model light curve calculations presented in 
this section the following parameters were chosen. The model 
atmosphere is divided into 128 layers. The radiative transfer is 
solved including atomic lines of the Kurucz atomic data line list. 
The number of wavelength points used for the solution of the ra- 
diative transfer is about 2400. For the inner boundary condition 
of the radiation, we use the nebular boundary condition, which 
means that the radiation is just passing through the inner empty 
region, therefore all the radiation is produced by the ejecta itself, 
there is no inner lightbulb. The outer boundary condition is zero. 

For the first LTE light curve calculation, the W7 deflagration 
explosion model is used to obtain model light curves of SNe la 
in different photometric bands. The first point of the light curve 
was calculated at three days after the explosion. The method to 
obtain the theoretical light curves is described in section |4~T1 

Before we present our resulting model light curves, we 
present plots of various terms that change the energy density for 
several moments in time. In Fig. [8] the energy change due to 
the expansion is shown for several moments in time. The energy 
change is decreasing for later days in the light evolution. The 
additional energy input due to the radioactive decay is shown in 
Fig. [9] With increasing time, the amount of energy input into the 
atmosphere decreases. It also shows, that the additional energy 
is located in the inner part of the envelope, where the 56 Ni has 
been produced. Fig. [TO] shows the energy loss due to the emis- 
sion of radiation. The temperature structures of the atmosphere 
for several moments in time are shown in Fig. Q~T] 

Figure[T2]shows the LTE SN la model light curve of the W7- 
based explosion model in the optical V band. The theoretical 
light curve accurately reproduces the observed light curves of 
two SNe la. The steep rise of the model light curve beginning at 
three days after explosion is in agreement with the observed light 
curves. The maximum of the W7 -based model light curve seems 



Fig. 8. The energy change due to the expansion for several days. 
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Fig. 9. The energy change due to the additional energy input for 
several days. 



to be later than that of the observed light curves. At 20 days after 
the explosion, the model light curve has its maximum, while the 
maximum of the observed light curves is around 17 days after the 
explosion. After maximum, the decline of the light curve of the 
W7-based model well reproduces the observed light curve. Even 
up to the later phase at 50 days after the explosion, where the 
atmosphere becomes significantly thinner, the fit to the observed 
light curves is quite accurate. 

The theoretical light curve in the ultraviolet U band is shown 
in Fig. [13] Only an observed light curve of SN 1999ee is avail- 
able. The observational data are scattered. The rise in the begin- 
ning as well as the maximum phase is well represented by the 
model light curve. However, the decline of the theoretical light 
curve seems to be too steep. This same effect is present in the 
model light curve of the B band, which is shown in Fig. [14] The 
first days of the model light curve are too bright compared to 
both observed SN la light curves. The maximum phase of the 
model light curve is in good agreement with the observed ones. 
At day 50, the model light curve becomes brighter than the ob- 
served light curves. 

In Fig. [15] a plot of the model light curve of the R band is 
shown. The steep rise in the beginning and the maximum phase 
of the observed SN la light curves is well represented by the 



6 



D. Jack et al.: Theoretical light curves of type la supernovae 



.5 -5.0x10 s 



-1 .0x10 1 



ET -1.5x10" 




1 x10 4 
expansion 



day 10 
day 20 
day 30 
day 40 
day 50 



2x10 4 
'elocity [km/s] 



3x10 4 



-19.5 



19.0 



-18.5 



light curve in V band 




20 30 
days since explosion 



Fig. 10. The energy change due to the radiation for several days. 
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Fig. 11. Temperature profiles as a function of time. 



Fig. 12. LTE model light curve of the W7 explosion model in 
the V band compared to two observed SN la light curves of SN 
1999ee and SN 2002bo. 



-19 



-17 



light curve in U band 
















/i D 
B 


□ SN 1999ee -| 






\ D - 






\ □ 






D □ 






\ □ 






\a a 






\ □ 












\B D n : 






\ □ □ 






\ n Id 






\ □ 






\ □ — 






\ § 






\ pBp D l 






X HnB □□ : 






\ £ 






\. □- 










10 20 30 40 50 

days since explosion 



computed model light curve. However, the theoretical light curve 
fit becomes worse for later phases. The luminosity of the model 
light curve seems to rise again at around day 25 after the explo- 
sion. Up to day 45, the model light curve has a second bump, 
which is not observed in the light curves of SN 1999ee and SN 
2002bo. In the infrared I band, the decline after the maximum 
phase is missing, as shown in Fig. [16] As in the R band, the rise 
in the beginning and maximum are well represented in the model 
light curve. However, at maximum, the luminosity of the SN la 
model light curve rises further, which is not seen in the observed 
light curves of SN 2002bo and 1999ee. Around day 30, the dif- 
ference between model and observed light curve in the I band are 
about 1 mag. Up to day 50, the model light curve declines, while 
the observed light curves show their second maximum around 
40 days after explosion. 

In Fig. [17] we compared the relation U-B of the theoretical 
and observed light curves. The same comparision for the B-V 
relation is shown in Fig. [18] 

4.3. Dynamical models 

We have three different results of explosion calculations of SN 
la events. The structure of these models are results from hydro- 
dynamical explosion calculations. One is the W7 deflagration 



Fig. 13. U band light curve of the LTE model light curve com- 
pared to SN 1999ee. 

model of No moto et al.l d 19841) . The W7 dynamical model as- 
sumes that the ongoing explosion is a deflagration. The flame 
propagates with a velocity lower than the speed of sound out- 
wards. Another possible explosion model is the delayed detona- 
tion model. The explosion starts with a deflagration which even- 
tually proceeds to a detonation. We used two different dynamical 
models named as DD25 and DD16. For more det ailed informa- 
tion ab out the used delayed detonation models see Hoflich et al. 
(12002b . 

The resulting optical model light curves of the V band are 
plotted in Fig. [19] The light curve of the DDI 6 is fainter than the 
other ones. This is due to the fact that in this explosion model 
less 56 Ni is produced. The W7 and DD25 model are quite similar. 
DD25 is slightly brighter than W7. The shape of the light curves 
are in good agreement with the observed ones. The maximum of 
the DDI 6 is later than the one of the W7 and DD25. 

In Fig.[20]the light curves in the I band of the three dynami- 
cal models are shown. The infrared light curves deviate strongly 
from the observations. As in the optical light curve, the DDI 6 
light curve is fainter than W7 and DD25. DDI 6 also has only 
one maximum and a decline. However, this maximum is later 
than the observed ones. The DD25 and W7 light curves are still 
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Fig. 14. LTE light curves of the W7 explosion model. The B Fig. 17. Comparison of U-B for the theoretical and observed 
band model light curve is too bright during the first few days. light curves. 
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Fig. 15. The theoretical light curve in the R band seems to rise Fig. 18. Comparison of B-V for the theoretical and observed 
after a first decline after the maximum phase. light curves. 
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Fig. 16. LTE light curves of the W7 explosion model. For the I Fig. 19. V band light curves for three different dynamical models 

band the model light curve has no distinctive maximum. At 30 W7, DDI 6 and DD25. The DD25 dynamical model is fainter 

days after the explosion the model light curve is way too bright because less nickel was produced during the explosion. The W7 

in comparison to the observed ones. and DD15 dynamical model light curves are quite similar. 
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Fig. 20. 1 band light curves for the dynamical models W7, DD 16 
and DD25 are compared to the two observed SN la light curves. 



rising after the observed light curves have reached their maxi- 
mum and are declining. The theoretical light curves in the in- 
frared must be improved for all explosion models. 

5. Light curves of NLTE models 

The LTE model light curves in the V Band and most other 
bands are in quite good agreement with observed SN la light 
curves. However, in the near-infrared of the I band, the model 
light curves need improvements especially for the later phase 
to fit the observational light curves more accurately. So far, the 
SN la model atmosphere is considered to be in LTE during the 
whole evolution time. In this section, we calculate the model 
light curves with the assumption of an atmosphere which is not 
in LTE. The computation of model light curves with SN la model 
atmospheres in NLTE requires a huge amount of computation 
time. At first, we model light curves with atmospheres in NLTE 
that are computed with LTE temperature structures. More real- 
istic NLTE model light curves with a temperature structure that 
adapts to NLTE conditions are computed to investigate the NLTE 
effects on the model light curves. 

The first approach to compute NLTE model light curves is 
to consider the atmosphere to be in NLTE, but use a fixed LTE 
temperature structure. For this computation of NLTE model light 
curves, we used the calculated radiative equilibrium LTE tem- 
perature structure of the W7 deflagration model. We keep this 
temperature structure constant and perform 20 iterations to let 
the NLTE converge, which is mainly the occupation numbers 
of the species that are considered for the NLTE calculation. We 
considered the following species for the calculations in NLTE: 
H I, He I, He II, C I-III, O I-III, Ne I, Na I, Mg 1-TR, Si I-III, 
S I-III, Ca II, Ti II, Fe I-III and Co II. These are the species that 
are most abundant in SN la atmospheres and mainly contribute 
to the spectrum. There is only little H and He in SN la, and in 
W7, there is none. However, considering H I, He I and He II 
in NLTE does not cost much computation time. The advantage 
of this approach to an NLTE model light curve is that no tem- 
perature iterations have to be performed. This reduces the com- 
putation time significantly, although about 200,000 wavelength 
points are calculated instead of 2,400 in case of LTE. 

In Fig. [2X1 the SN la model light curve in the V band of an 
NLTE calculation with an LTE temperature structure is shown. 
The NLTE light curve is fainter than the LTE light curve. During 



Fig. 21. Model light curves of the W7 explosion model in the 
V band. The NLTE model atmosphere has an LTE temperature 
structure. 
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Fig. 22. Model light curves of the W7 model in the U band. The 
NLTE model atmosphere has an LTE temperature structure. 



the maximum phase there is about 0.4 mag difference between 
both light curves. We have to point out that for our NLTE 
models, the energy is not conserved because the radiation does 
not thermalize within the envelope. After maximum the NLTE 
model light curve approaches the LTE light curve and agrees 
with the observed light curves as well as the LTE light curve. 
However, the NLTE model light curve in the V band does not 
agree with the observed light curves very accurately. With the 
assumption of an atmosphere in LTE, a better fit is obtained, al- 
though NLTE is more accurate. 

The NLTE model light curve of the U band is shown in Fig. 
1221 The NLTE light curve shows almost no deviations from the 
LTE light curve. During the later phase the steep decline is also 
present in the NLTE model light curves. In Fig. [23] the NLTE 
model light curve in the B band is presented. The NLTE model 
light curve is slightly fainter than the LTE light curve. The shape 
of the NLTE light curve seems to be the same as for the LTE 
light curve. The NLTE model light curve is also an accurate fit 
to the observed light curves. 

In Fig. [24] the NLTE model light curve of the R band is 
shown. The luminosity of the NLTE light curve is fainter than 
the LTE light curve. It has almost the same shape. But for the 
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Fig. 23. Model light curves of the W7 model in the B band. The Fig. 26. NLTE and LTE model light curves of the W7 explosion 
NLTE model atmosphere has an LTE temperature structure. model in the V band. 
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Fig. 24. NLTE and LTE model light curves in the R band of the 
W7 model. The NLTE model atmosphere has an LTE tempera- 
ture structure. 
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Fig. 25. NLTE and LTE model light curves of the W7 model. The 
NLTE model atmosphere has an LTE temperature structure. In 
the I band the NLTE light curve is somewhat of an improvement. 



phase after 25 days after the explosion, the NLTE model light 
curve rises again. The assumption of NLTE does not improve 
the fit to the observed light curves of SN 1999ee and SN 2002bo. 
The infrared NLTE light curve for the I band is shown in Fig. |25] 
During the maximum phase, the NLTE light curve is fainter than 
the LTE light curve. A distinctive maximum is missing in the 
NLTE model light curve. For the phase between day 15 and day 
25, the NLTE model light curve fits the observed SN la light 
curves very accurately. Here, there is a significant improvement 
compared to the light curve of the LTE model atmosphere. But, 
after day 25 the NLTE model light curve starts to rise and be- 
comes too bright. This is the same problem that already emerged 
in the LTE light curve. 

5.1. NLTE atmosphere structures 

We perform a more realistic NLTE calculation of the SN la 
model atmosphere evolution. The temperature structure is now 
changing and adapting to the new conditions of NLTE. The cal- 
culation of an NLTE model light curve takes a huge amount 
of computation time. Considerably more wavelength points are 
needed for the calculation of the solution of the radiative trans- 
fer. For all the species considered in NLTE, the rate equations 
have to be solved. Note that a time step in the NLTE calcula- 
tion is not a real time step. The rate equation changes the energy 
of the atmosphere, but this is not included in the energy solver. 
However, it is adequate as the goal is to obtain a temperature 
structure, where the atmosphere is in radiative equilibrium. 

In a first try to calculate a more realistic NLTE light curve, 
numerous species up to calcium are considered to be in NLTE. 
These are the species H I, He I, He II, C I-III, O I-III, Ne I, 
Na I, Mg I-III, Si I-III, S I-III and Ca II. Higher species are 
neglected because they have more levels, which would increase 
the computation time significantly. Nevertheless, the computa- 
tion of an NLTE light curve needs significantly more time than 
LTE. For this computation of NLTE model light curves, about 
200,000 wavelength points need to be calculated, compared to 
about 2,400 wavelength points for a pure LTE light curve cal- 
culation. To obtain an NLTE model light curve in a reasonable 
time, we started the calculation at day 10 after the explosion 
and used the LTE temperature structure as the initial atmosphere 
structure. The main focus is to check if the infrared light curves 
during the later phase can be further improved. 
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Fig. 27. NLTE and LTE model light curves of the W7 explosion 
model in the I band. 
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Fig. 28. Light curves in the I band. The light curve with variable 
scattering is fainter than the LTE light curve. 



In Fig. the NLTE and LTE light curves of the W7 defla- 
gration model in the V band are shown. The maximum phase is 
not well fitted by the NLTE model light curve. The LTE light 
curve fits the observed light curves better. At day 20 the NLTE 
and LTE light curves are almost the same. Compared to the 
NLTE light curve obtained with the LTE temperature structure, 
there are only small differences to the NLTE calculation where 
the temperature structure adapts to the NLTE conditions. 

In Fig.|27l the I band model light curves of NLTE and LTE 
are shown. During the maximum phase, the NLTE light curve is 
fainter than the LTE light curve. Between day 15 and day 25 the 
NLTE model light curve fits the observed light curves quite ac- 
curately. Here, the use of NLTE improves the fit to the observed 
light curves. However, we obtained this improvement also with 
the NLTE model light curve calculated with the LTE temperature 
structure. At day 25 the NLTE model light curve starts to rise 
again. Although the consideration of NLTE improves the model 
light curve in the I band, the problem with a rise in brightness 
after the maximum remains. 



5.2. Line scattering 

The fits to the infrared light curves need to be improved, as they 
are too bright during the phase of 20 to 50 days after the explo- 
sion in all models. We found that a change in the ratio of line 
scattering to absorption in the LTE case improves the fits to the 
I band light curve. The resulting light curve is shown in Fig. 
1281 The line scattering is here roughly inversely proportional to 
the density and was for the figure "hand tuned" to produce a 
good match. As the density becomes lower due to the ongoing 
expansion, the scattering becomes more important, as expected. 
A more detailed investigation of the IR line scattering and its 
effects on SN la light curve modeling will be a topic of future 
work. To treat this effect more accurately, very detailed NLTE 
model atoms are required that include the important infrared 
lines. In the atomic data we have presently available, these lines 
are mostly "predicted" lines and were not included in the our 
current NLTE model atoms. This result indicates that an accu- 
rate treatment of line scattering and NLTE effects in the IR lines 
is very important to model and interpret the IR light curves of 
SNe. 



6. Conclusion 

We have presented a new approach for a energy solver which 
can be applied to expanding SN la envelopes to compute SN 
la model light curves. Test calculations confirm that the imple- 
mented code works properly. We applied the energy solver to 
calculate SN la model light curves during the free expansion 
phase. 

We present the first PHOENIX light curves of type la super- 
novae in different photometric bands. At first, we solved the ra- 
diative transfer equation with the assumption of the atmosphere 
being in LTE. The model light curves were in decent agreement 
with the observed ones at early times. However, in the infrared 
I band, the theoretical l i ght cu rves do not fit the observed light 
curves. iBlinnikov et all (2006) had the same deviations in their 
model light curves compared to observations. Their light curve 
in the I band is also too bright and furt her rising, e ven after the 
maximum in observed light curves. iKasenl (120061) presented a 
detailed study of SN la light curves in the near IR. His I band 
light curve fits the observations better. He also showed how vari- 
able and dependent on different parameters the light curves in 
the near infrared can be. 

Three different dynamical models were compared to ob- 
served SN la light curves. The delayed detonation model DD 
16 is unlikely to be the best explosion model because the model 
light curves are too faint to fit the observed light curves due to 
the low 56 Ni mass. The results of the W7 model are the best fit 
to the observed light curves. 

For the later phase of the light curves, especially in the in- 
frared, we need the atmosphere to be considered in NLTE. We 
showed that NLTE model light curves computed with an LTE 
temperature structure leads to improvements in the model light 
curves. For future work more investigation of the infrared model 
light curves and the NLTE needs to be performed. The focus 
has to lie on the influence of line scattering on the shape of the 
light curves. It would also be interesting to look in detail on the 
spectral evolution during the free expansion phase. Further im- 
provements may be achieved with a full 3D radiative transfer 
and energy solver. 
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